Hybridization of spin and plasma waves in Josephson tunnel junctions with a 

ferromagnetic layer. 
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We study dynamics of tunnel Josephson junctions with a thin ferromagnetic layer F 
[superconductor- insulator- ferromagnet-superconductor (SIFS) junctions]. On the basis of derived 
equations relating the superconducting phase and magnetic moment to each other we analyze col- 
lective excitations in the system and find a new mode which is a hybrid of plasma-like and spin 
waves. The latter are coupled together in a broad range of parameters characterizing the system. 
Using the solution describing the collective modes we demonstrate that besides the Fiske steps 
new peaks appear on the FV characteristics due to oscillations of the magnetic moment M in the 
ferromagnetic layer. Thus, by measuring the FV curve of the SIFS junctions, one can extract an 
information about the spectrum of spin excitations in the ferromagnet F. 

PACS numbers: 74.50.-l-r, 03.65.Yz, 74.20.Rp, 85.25.Cp 



Discovery of the Josephson effects in 1962 was an im- 
portant step in the development of the condensed matter 
physics [H. Josephson has predicted that the dc cur- 
rent of Cooper pairs with a magnitude of order of the 
qusiparticle current can flow in the absence of the volt- 
age through a thin insulating layer separating two super- 
conductors S. In the presence of a voltage V the super- 
current oscillates in time with the Josephson frequency 
Ljj = 2eV/h. 

On the basis of equations describing electrodynamics of 
the Josephson tunnel junction [superconductor-insulator- 
superconductor (SIS)] many fascinating phenomena were 
predicted [li0]. For example, small perturbations of the 
phase difference may propagate in SIS junctions in a 
form of Josephson "plasma" waves with the spectrum 
uP' = Oj(l -I- {MjY). Large perturbations exist in a form 
of solitons carrying the magnetic flux quantum (fluxons 
or antifluxons). The interaction of the "plasma" waves 
with oscillating Josephson currents leads to resonances 
and peculiarities on the I-V characteristics (Fiske steps) 
i- 

The Josephson junctions (JJ) based on conventional 
superconductors are widely used in practice as genera- 
tors, the most sensitive detectors of magnetic fields and 
ac radiation in a wide range of the frequency specr- 
tum, etc [l]. Recently, a considerable progress has been 
achieved in practical applications of the JJs based on high 
Tc superconductors [31]. 

Replacing the insulating (I) layer by a ferromagnet (F) 
one obtains SFS junctions that have been under intensive 
study during the last decade [1, |3, 0, 01 ■ interest- 
ing effects like the so-called tt— state (negative critical 
Josephson current jc) or new type of superconducting 
correlations (odd triplet superconductivity) may arise in 
JJs of this type [1, @j S| • Although the main attention 



was paid to the study of the dc Josephson effect, some ac 
effects in SFS or SmS (m means a magnetic nanoparticle) 
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junctions have been considered. 

At the same time, all these studies deal with the SFS 
junctions where effects like those in tunnel JJs are ab- 
sent. Experimental studies of the SIFS junctions with an 
insulating layer I have started quite recently (see [l3| and 
references therein). In these tunnel JJs one can observe 
collective modes and other dynamical phenomena. Note 
that the effect of superconductivity on the ferromagnetic 
resonance has been studied experimentally on S /F bilay- 
ers [12]. We have no doubts that proper experiments will 
be performed in the nearest future and theoretical de- 
scription of dynamical effects in SIFS or SFIFS tunnel 
j mictions is in great demand. 

By now, no theory for dynamical effects in SIFS or 
SFIFS tunnel junctions has been suggested, although one 
can expect really new effects. One of the basic proper- 
ties of the JJs is that they are very sensitive to even very 
small variations of the magnetic field. So, one can eas- 
ily imagine that dynamics of the spin degrees of freedom 
leading to variations of the magnetic field should essen- 
tially affect electrodynamics of tunnel SIFS junctions. 

In this paper, we develope theory describing dynamical 
effects in SIFS or SFIFS junctions. We demonstrate that 
dynamics of the phase difference and of the magneti- 
zation M is described by two coupled equations which 
determine dynamics of both the spin and orbital degrees 
of freedom. One of these equations governs the spatial 
and temporal evolution of the phase (/?, while the other 
describes the precession of the magnetization M. Using 
these equations we investigate collective modes in Joseph- 
son junctions with F layers and show that these modes 
are coupled spin and plasma-like Josephson waves. 

Although the influence of the superconducting con- 
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densate on spin waves in magnetic superconductors was 
studied by Braude and Sonin Q some time ago, the 
tunnel JJs with a ferromagnetic layer is a more interest- 
ing system because besides the spin waves another mode 
(plasma-like Josephson waves) exists in these systems. It 
is our main finding that in such systems the both types 
of the excitation hybridize forming a new collective mode 
consisting of coupled oscillations of the magnetization 
and the Josephson current. Moreover, the spin waves can 
be excited and recorded simply by passing a dc current 
j through the junctions, and measuring the I-V char- 
acteristics of the junction. We demonstrate below that 
peculiarities on the I-V curve (Fiske steps) bear infor- 
mation about the spectrum of the system, which can be 
considered as the basis of a new type of spectroscopy of 
spin excitations in a ferromagnet. 

The hybridization of the charge and spin excitations 
occurs because the contribution to the magnetic moment 
in the system comes from both orbital electron motion 
and exchange field. The orbital currents are affected by 
spins of the ferromagnets and vice versa. The coupling 
between the plasma-like and spin waves in SIFS junctions 
may show up in additional resonances and peaks on the 
I-V characteristics. 

Having discussed the physics of the SIFS junctions 
qualitatively, let us present now the quantitative theory. 
We consider for simplicity an SFIFS junction (the ob- 
tained results are valid also for an SIFS junction shown 
schematically in Fig. la) and calculate the in-plane cur- 
rent density j± as a function of the magnetic induction 
B± in the plane parallel to the interface. This cur- 
rent can be expressed through the vector potential Aj_ 
(V X A = B) and gradient of the phase V^x 
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where is the London penetration length and <f>o — 
hc/2e is the magnetic flux quantum. 

Writing Eq. Ilj we imply a local relationship between 
the current density j_L and the gauge invariant quantity 
in the brackets, which is legitimate in the limit fcA^ ^ 1, 
where k is the modulus of the in-plane wave vector of per- 
turbations. Subtracting the expressions for the current 
density, Eq. written for the right (left) supercon- 

ductors from each other we find a change of the current 
density [j^] across the junction 




FIG. 1: A) Schematic picture of the system under consider- 
ation. The thick line between the S and F layers means an 
insulating layer; B) Spectrum of coupled spin and plasma-like 
modes. On the vertical and horizontal axes we plot the quan- 
tity (lo/Qj)'^, and (klj)^, respectively. The following param- 
eters are chosen: (Q.j/flM)'^ — 2, (Ij/Im) ~ 1 and s = 0.05 



We proceed by solving the London equation in the su- 
perconductors with boundary conditions determined by 
the change of the current density (see, e.g., [l^l)- A so- 
lution of this equation for S films with the thickness ex- 
ceeding Xl takes the form 



B^(z) = { 



$0 „ 27rdF-, , . 
nzXVj_ip — M^}exp(- 
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(3) 

Eq. ^ describes the penetration of the magnetic induc- 
tion into the superconductors (|z| > dj?). It is well known 
that the magnetic field penetrates the J J in a form of fiux- 
ons P, Q • Eq. ([3]) supports this picture for the SFIFS 
junctions. Integrating Eq. ([3]) over z and x and adding 
the magnetic flux in the F film — AndpLxMy (the 
magnetic moment M±^ is assumed to be oriented in the 
y direction), we come to flux quantization: $ — [(^^n)^ 
where n = [(^]a;/27r is the number of fluxons and [lp\x is 
the phase variation on the length Lx^ which is supposed 
to be an integer multiple of 27r. 

In order to obtain an equation for the phase difference 
(/?, we use the Maxwell equation (V x B)^ — {An/c)jz in 
the superconductors and the standard expression for the 
Josephson current. Then, we obtain 



ic/47T)Xl^{47Tdp{n,xM^) + ($o/2^)Vi^} (2) 



where dp = dp (or 2dp) in the case of a SIFS (or SFIFS) 
junction, dp is the thickness of the F film which is as- 
sumed to be smaller than the London penetration length 
Xl- This assumption allows one to neglect the change of 
A^ caused by Meissner currents in the F layer and to 
write the change of the vector potential A^ in the form 
[A_l] = dp{nzXB±_) with B_l = 4nM±_ + Hj^. 
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(4) 

where flj — {2ejc/C\jhy/^ is the Josephson "plasma" 
frequency, 7 = (i?nCn)^^, Cn = e/Aird and i?n are the 
capacitance and resistance of the junction per unit area, 
d is the thickness of the insulating layer, /j = u}/r2j, 
vj = cy^ d/2eXL is the velocity of the plasma wave prop- 
agation (Swihart waves). The first term on the r.h.s. of 
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Eq. (|4]) , r] = j/jc, is the normalized current through 
the junction. A simpler equation for the stationary case 
has been written previously in Ref. [l5| . where a similar 
system with a multidomain ferromagnet was considered. 

The dynamics of the magnetization in the F layer is 
described by the well known equation with account for 
the magnetic induction due to the Meissner currents (see 
for example Ref. [3]) 



9Mj 
dt 



-4T:ga{l-lli\7l){MxM^)+gMxB^ (5) 



where g is the gyromagnetic factor, a is a parameter 
related to the anisotropy constant /? = (a — 1), is a 
characteristic length related to the spin waves. 

We further substitute 'B±p = 4:n'M± + B±s{dF) into 
Eq. (O, where Bj^^ is the magnetic induction in the F 
layer and 'B±g = B^(z dp) is given by Eq.®. Finally 
we come to the equation 
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dt Mo 

(6) 

where — '^T^gM^jS — A:-KgMQ{a — 1) is the resonance 
frequency of the magnetic moment precession {a > 1), 
s — dp/{2f3XL),l\i = (a/{a — \))l\j. One could also take 
into account a damping adding to the r.h.s. of Eq. ^ 
the term 7a/(M x dl>A./dt)/M {-fM is the dimensionless 
Gilbert constant). 

Eqs. ([H ini) are the final equations fully describing 
dynamics of the SFIFS junctions. The most interest- 
ing effect that we will obtain now is the hybridization of 
charge and spin excitations resulting in a new collective 
hybridized mode. 

a) Hybridized mode. Let us consider the simplest mon- 
odomain case when the magnetization Mq is normal to 
the interface, so that in equilibrium Bq — 0. Small per- 
turbations near the equilibrium result in a precession of 
the magnetic moment M and in a variation of the phase 
difference ip in time and space. Representing M for small 
deviations mj^ as M —Mon2+m.j_, we can linearize Eqs. 
(|6]) and ([4]). A nonzero solution of these linearized equa- 
tions exists provided the determinant of these two equa- 
tions equals zero. Writing the perturbations in the form 
of plane waves (m^ ip exp{iujt + ikrj^)) and set- 
ting the determinant to zero we come to the dispersion 
relation 



(7) 



where illj{k) = ^m{1 + s 

{kljY). For simplicity, we neglected the damping setting 
7m 0. 

Eq. ([7|) describes the spectrum of the hybridized mode. 
This mode decouples into the spin and charge excita- 
tions only in the limit when the right hand side can be 




FIG. 2: Correction to the I-V characteristics due to interac- 
tion of Josephson oscillations and collective modes. The nor- 
malized voltage is defined as Vnorm ~ i^v/^j- Solid curves 
correspond to a small parameter s = 2(dF / P^l) = 0.05 
(negligible effect of the F layer), point curves correspond to 
s = 0.4. The plots are presented for the following parameters: 
A) Q.m/^j = 2, n„lj = 2; B) fiAz/fij = 1, t^nlj = 4. The 
damping coefficients are '^j/Q.j — 0.2, jm /^j ~ 0.1. 



neglected. In this case the spin waves with the spec- 
trum flM^k) and the plasma-like Josephson waves with 
the spectrum fij(fc) exist separately. In the general case 
we have the coupled spin wave and plasma-like modes in 
the system. The most interesting behavior corresponds 
to the case Hm < ^j- In this situation, the two branches 
of the spectrum cross each other in the absence of the 
coupling and the coupling leads to a mutual repulsion of 
these branches. In Fig. lb we show the dependence u^k) 
just for this case. Note that the frequency of the con- 
sidered collective modes remains unchanged by inversion 
of the magnetization direction, M <-> —M, that is, in a 
multidomain case we would obtain the same dependence 
uj{k). 

The collective modes in the Josephson junctions can 
be excited by the internal Josephson oscillations. It is 
well known that the interaction of Josephson oscillations 
with plasma- like waves in a tunnel Josephson junction in 
the presence of a weak external magnetic field leads to 
Fiske steps on the TV characteristics (see, e.g., Refs. [2], 
,16.] )■ In the next paragraph we consider the modification 
of Fiske steps in tunnel SFIFS or SIFS junctions due to 
the hybridization of the collective modes using Eqs. ([H 

ED. 

b) Fiske steps. We assume again that the magnetiza- 
tion M is oriented perpendicular to the interfaces and a 
small external magnetic field H^. is applied parallel to the 
films. It is important to have in mind that the magneti- 
zation Mq is of the order of lOOGauss or larger, whereas 
He is about a few Gauss. Therefore one can neglect 
the in-plane magnetization My = i7e/47r/3 in comparison 
with Mq. We assume also that the length of the junction 
along the x— axis is shorter than the Josephson length Ij. 
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In this limit the phase can be represented in the form 
(see Refs. [l^): = ipo{x,t) + ■>p{x,t) with ipo{x,t) = 
flyt + khx and fly = 2eV/h, — 47rAi7?e/^o- Lin- 
earizing Eq. ^ we obtain the equation for ip 



,2^ 

' dx^ 



- sm (fo ■ 



cdp drUy 
2\Ljc dx 



(8) 



The dc current is 77 = {2eV^/h)/Q,j + 
(^(x, i) cos (/3o(a;, ^)), where the angular brackets de- 
note the averaging in space and time. The derivative 
dniy/dx can be found from Eq. ([6]) with ip replaced 
by "0. The first term on the right-hand side in Eq. ([8|) 
plays a role of an external force oscillating in space and 
time and acting on a resonance system. This equation 
should be solved using boundary conditions dtp/dx = 
at a; = ±Lx/2. Carrying out these calculations we obtain 
a contribution Sr] = t) cos (/3o(a:, t)) to the current 

originating from the collective modes 



e^fj cos(26ly) - cos(26l^). 

0^ {el^el)sm{2ev) 



} (9) 



where D ^ {e'y - eij)a'y / O'y , a'y = (f^f, _ tjjny)/nj, 

9h = KhLx/2,9v = KvLx/2,Ky = ay/ly, ly = + 

snlJi^l - = -l- s). For simplicity 

we consider the limit Ij » Im- 

In Fig. 2 we plot the dependence 5rj{V) vs the normal- 
ized voltage Jly/Oj for different values of the parameter 
s = 2'KdF / (PXl)- In order to avoid a divergence, we 
took into account a finite damping in the F layer replac- 
ing (f]^ - nlj^) by {ny{ny - i-iu) ~ ^Its)- 

In the limit of small values of s (s — > 0) this depen- 
dence describes Fiske steps in the conventional Josephson 
SIS junction. Solid lines in Figs. 2 correspond to small 
values of s (s = 0.05) when the influence of the F layer 
is negligible and point lines correspond to larger values 
of s (s = 0.4). One can see from Fig. 2a that near the 
frequency of the magnetic resonance, ily « ^ms = SOj 
there is a double peak (point curve), which is absent in 
junctions with a small s (solid curve). Fig. 2b shows 
that the presence of the F layer leads not only to a shift 
of maxima in the dependence 6rj{V) but also to a change 
of the overall form of this dependence. 

To conclude, we have developed theory describing elec- 
trodynamics of the SFIFS or SIFS tunnel Josephson junc- 
tions. In the framework of our approach dynamics of 
these systems is fully described by Eqs. ([310. Solving 
these equations we found new effects related to hybridiza- 
tion of charge and spin degrees of freedom. We demon- 
strated that the different branches - Josephson plasma- 
like and spin wave branch - of the spectrum repel each 
other. The hybridization of the collective modes results 
in new interesting peculiarities of the I — V character- 
istics. We have demonstrated that the positions, shape 



and number of the Fiske steps arising in the presence of 
a weak external in-plane magnetic field and an applied 
voltage V change. In particular, new peaks appear due 
to excitation of magnetic modes. We believe that on the 
basis of the SIFS junctions one may construct a new type 
of spectrometers that would allow one to extract an in- 
formation about spin excitations in thin ferromagnetic 
layers. It seems that an experimental observation of the 
effects predicted here is not very difficult and can be per- 
formed by measuring I — V characteristics, spectra of 
collective excitations, etc., using standard experimental 
methods developed for studying Josephson junctions. 
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